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We report on the long-term evolution of the spin period of the symbiotic X-ray pulsar GX 1+4 and 
a possible interpretation within a model of quasi-spherical accretion. New period measurements 
from BeppoSAX/WFC, INTEGRAL/ISGRI and Fermi/GBM observations have been combined 
with previously published data from four decades of observations. During the 1970's GX 1+4 was 
spinning up with the fastest rate among the known X-ray pulsars at the time. In the mid 1980's it 
underwent a change during a period of low X-ray flux and started to spin down with a rate similar 
in magnitude to the previous spin up rate. The spin period has changed from ^HOsto^^Os 
within the last three decades. Our results demonstrate that the overall spin down trend continues 
and is stronger than ever. We compare the observations with predictions from a model assuming 
quasi-spherical accretion from the slow wind of the M giant companion. 
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1. Introduction 

Accreting X-ray pulsars are highly magnetized neutron stars in a binary system accreting mat- 
ter from their companion star. The mass transfer can take place via Roche-Lobe overflow, strong 
stellar winds for giant stars or the Be emission mechanism. These accreting pulsars radiate pre- 
dominantly in the X-ray band, and the radiation is modulated by the stellar rotation of the pulsar. 
For a review on this subject, see e.g., fli~5j]. 

GX 1+4 is an accreting X-ray pulsar discovered in 1970 by a balloon X-ray observation at 



energies above 15 keV showing pulsations with a period of about two minutes ([|13j]). It was one 
of the brightest X-ray sources in the Galactic center region. The composite emission spectrum of 
GX 1+4 indicated that the object was almost certainly a binary system, consisting of a symbiotic 
red giant and a much hotter source ([[|]). [0 confirmed the optical companion to be the M giant 
V2116 Oph. 

More recently, observations have shown a well-determined, 1161 day period, single-line spec- 



troscopic binary orbit ([12]). The inclination of the orbit is, however, unknown. Therefore, on the 
basis of similar characteristics seen in other X-ray binary systems the mass of the neutron star is 
inferred to be ~ 1.35 M®. The mass function of the system, combined with the abundances of the 
M giant star and the assumed mass for the neutron star, indicate a mass of ~ 1 .2 M & for the M giant 
star. This implies that the M giant star is a first ascent giant which does not fill its Roche Lobe. 
In consequence, the neutron star is capturing the slow stellar wind of its companion, which makes 



V2116 Oph quite different from other low-mass X-ray binaries (Q12[]). 

Currently there is no accurate estimation of the magnetic field of GX 1+4. Assuming the 
standard accretion disk theory, the magnetic field has been estimated to be 5 > 3.7 x 10 13 G([|3|]), 
which is among the largest measured for accreting X-ray pulsars. However, there have been de- 
tected some signatures of the existence of CRSFs in the X-ray spectra of the neutron star which 
points to a value of B ~ 10 12 G for the magnetic field ([[7|]). 

Therein, GX 1+4, is a unique pulsating X-ray binary in many aspects. Among its singularities, 
it is the first low mass X-ray binary with an M giant optical companion, the first symbiotic X-ray 
pulsar ever observed, could have the highest magnetic field of any known accreting X-ray pulsar 
and has a peculiar spin behavior which explanation is the aim of this work. 

2. Results 

The main result obtained during this work is the plot shown in Fig. [T] which currently is the 
most updated and complete series of pulse periods of GX 1+4. The available measurements of 
the pulse period of the neutron star of GX 1+4 span a time period of about 40 years giving us a 
unique insight in the pulse period evolution of this symbiotic X-ray pulsar. To obtain this result we 
have combined all measurements that we are aware of from the literature with the pulse periods 
determined from BeppoSAXAVFC, INTEGRAL/ISGRI and Fermi/GBM data. 

The system has a peculiar spin history. During the 1970's GX 1+4 was spinning up with the 



fastest rate (P ~ 7.23 x 10 s/s) among the known X-ray pulsars at the time (see e.g., [15 ]). During 
that period the source had an X-ray luminosity between 9 x 10 38 and 2.25 x 10 38 erg s _1 ([p^]) 
(the luminosity value has been corrected to the distance range (3-6 kpc) ([0]) ). However, EXOSAT 
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Figure 1: Long-term pulse period of GX 1+4 ( [|1 OQ and references therein) spanning a time of about 40 
years. Note that the period increases from top to bottom. 



observations in 1983 revealed an extended low state with an X-ray luminosity of < 4 x 10 36 erg s _1 
(JTT|]). Around the same time GX 1+4 started to spin down with a rate similar in magnitude to 
the previous spin-up rate (see Fig. [I]) and continues spinning down up to date. Indeed, the spin- 
down rate measured by Fermi/GBM is the strongest spin-down rate observed to date (P ~ 1 .08 x 
10~ 7 s/s), and the pulse period of the pulsar has increased by about 50% during the last ~ 30 years 
to the largest value ever known for this source (~ 160 s) ([10]). 



Apart from the long-term plot, looking at the behavior along shorter time scales (Fig. ||), it is 
possible to observe that even when the X-ray flux has variations, pulse period evolution seems to 
be almost linear. 

The unusual long-term spin behavior of GX 1+4 has attracted considerable interest for many 
years. Studying the pulse period evolution in an accreting X-ray pulsar and relating it to, e.g., 
the luminosity changes, allows to test models of accretion and to gain insights on the physical 



processes taking place in this kind of systems. In [|10J we extend the investigation of the spin period 
history of GX 1+4 with new observations obtained by BeppoSAX/WFC, INTEGRAL/IS GRI and 
Fermi/GBM, we correlate X-ray flux variations with pulse period variations of the system using 
the same observations as well as those obtained with Swift/BAT and CGRO/BATSE, and provide a 
likely explanation for the spin history seen. 



3. Accretion models 

Standard accretion theory by [Q| assumes the formation of a prograde disk around the neu- 
tron star. Therefore, spin down of the pulsar is only possible in the quasi-equilibrium state, which 
implies the unusual magnetic field of B ~ 10 13 G. This prograde disk accelerates the neutron star 
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Figure 2: GX 1+4 data. Top left plot: Top: CGRO/BATSE daily lightcurve in the energy range 20-40 keV. Middle: Pulse 
periods derived from CGRO/BATSE data. Bottom: Residuals of the periods from a linear fit. Top right plot: Top: CGRO/BATSE daily 
lightcurve in the energy range 20-40 keV. Middle: pulse periods derived from BeppoSAX/WFC data. Bottom: Residuals of the periods 
from a linear fit. Bottom left plot: INTEGRAL/ISGRI data. Top: Average flux per pointing in the energy range 18 to 40 keV. Middle: 
20-40 keV pulse period. Bottom: Residuals of the periods from a linear fit. Bottom right plot: Top: Swift/BAT daily averaged light 
curve in the energy range 15-50 keV. Middle: Pulse periods derived from Fermi/GBM data. Bottom: Residuals of the periods from a 
linear fit. Note that in all these plots the period increases from top to bottom. 



through the angular momentum of the accreted matter. Thus, there is a positive correlation be- 
tween accretion rate (M) and the acceleration of the neutron star. Assuming M oc Fx, the positive 

fx In 

correlation v °< F x ([||]) should be observed in the data even during spin down. 

Transient disks with an alternating sense of rotation are known to form in numerical simu- 
lations in binary systems fed from stellar wind ([||]), hence, it is also possible the formation of a 
transient retrograde disk in GX 1+4 as it is a wind-fed source. Indeed, a retrograde accretion disk 
was proposed for the first time by [14] to explain the spin down observed in the neutron star of 
GX 1+4. This model assumes a retrograde disk around the neutron star decelerating the pulsar 
through the angular momentum of accreted matter, therein, the negative correlation -v« F^ 1 
should be observed in the data. Within this scenario the pulsar would necessarily be far from its 
equilibrium period. Therefore, a magnetic field of B ~ 10 12 G, instead of the high magnetic field 
predicted by standard accretion disk (B ~ 10 13 G), is required. 

In [10] we propose, for the first time, a third model to explain the behavior of GX 1+4, the 
quasi-spherical accretion model. Instead of a disk, this model predicts a quasi-static atmosphere 
around the neutron star which carries away angular momentum from the neutron star magneto- 
sphere but allows accretion (in contrast to the so-called subsonic propeller regime ([||])). As in the 
retrograde disk accretion model, this quasi-spherical scenario also predicts a negative correlation 
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Figure 3: The discrete correlation function between V obtained from Fermi/GBM GX 1+4 data and the 
Swift/BAT 15 - 50 keV X-ray flux shown in bottom right panel of Fig. 2. The minimum near zero lag 
implies a negative correlation between the pulse frequency derivative and the X-ray flux (||To||). 

between accretion rate and the acceleration of the neutron star during spin down, but within this 
scenario this correlation is weaker: — v °< F^ 1 ([jl8|]; [|l7|]). The formation of this quasi-static atmo- 
sphere is only possible for wind-fed sources like GX 1+4. This model does not require a magnetic 
field of B ~ 10 13 G as the pulsar does not need to be near its equilibrium period, but does not discard 
this unusual high magnetic field. 

4. Discussion 



The correlations between X-ray flux and pulse period evolution predicted by the models pro- 
posed in Section || can be tested by comparison with a regular sequence of X-ray fluxes and pulse 
periods which can be found in CGRO/BATSE and Fermi/GBM data (see Fig. |) for GX 1+4. 

A cross-correlation of v obtained from Fermi/GBM data and the contemporaneous X-ray flux 
measured by Swift/BAT finds a minimum near zero lag (see Fig. |J) implying a negative correlation 
between this two parameters for GX 1+4. Furthermore, this negative correlation has already been 
found previously in CGRO/BATSE data (i.e. [jl]]). As the standard accretion disk model predicts 
a positive correlation in every case, this negative correlation discards the standard accretion disk 
model as a possible interpretation of GX 1+4 behavior. 

Due to several reasons, including this negative correlation, different authors have tried to ex- 
plain GX 1+4 pulse period behavior assuming the formation of a retrograde disk around the neutron 



star by matter captured from stellar wind of the secondary M-giant star in GX 1+4 (e.g. []14[|, [gj, 
[0], [0], [0])- This model predicts a strong negative correlation between mass accretion rate and 
the acceleration of the pulsar (— V oc F^ 1 ). The correlation found between Fermi/GBM pulse peri- 
ods and Swift/BAT X-ray flux is — v °< F£ 03Q ([[K]]). Therein, the correlation found has the same 
sign as predicted by retrograde disk accretion model, but is about 3 times weaker than predicted. 
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Indeed, the negative correlation found in CGRO/BATSE data (— V °< F£ 0AS (||T|)) is also weaker 
than predicted by this model. Apart from this, the weak points of the retrograde disk interpretation, 
include the puzzling long-term (~ 30 years) stability of such a retrograde disk. 

This discrepancies between predicted and observed correlations and the puzzling long-term 
stability of the retrograde disk interpretation has led us to a new model, the quasi-spherical accre- 

3 11 

tion model. Whithin this model, there is no such a difference between predicted (— V °< F x ) and 



observed (— V °< F£ ([10]), —V °< i^ - 48 ([|IJ])) negative correlation, and it is possible to explain 
the observed correlation without the formation and maintenance, during such a long time period 
(~ 30 years), of a retrograde accretion disk ([10]). 



5. Conclusions 

New measurements by INTEGRAL/IS GRI and Fermi/GBM confirm the overall spin down 
of the neutron star in GX 1+4 observed over the last decades. This spin down is stronger than 
ever observed before with an increasement of the spin period of ~ 50% during the last ~ 30 years 



reaching the current value of ~ 160 s ([10]). The analysis of Fermi/GBM data reveals a nega- 
tive correlation between spin frequency derivative and the X-ray flux simultaneously measured by 
Swift/BAT monitor ([|10|]), which confirms the correlation of instantaneous spin-down torque with 
X-ray flux discovered by CGRO/BATSE ([[I]]) discarding clearly the standard accretion disk model 
for this system. We have shown that in GX 1+4 not only the retrograde disk accretion is a possible 
explanation but the quasi-spherical accretion onto the neutron star from the stellar wind of the M 



giant companion is likely to take place ([ 1C]). 
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